Collisional Properties of p-Wave Feshbach Molecules 
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We detected the formation of p-wave Feshbach molecules for all three combinations of the two 
lowest atomic spin states of 6 Li. By creating a pure molecular sample in an optical trap, we measured 
the inelastic collision rates of p-wave molecules. The elastic collision rate was measured from the 
thermalization rate of a breathing mode which was excited spontaneously upon molecular formation. 



Ultracold atomic gases with tunable interaction via 
Feshbach resonance have opened up possibilities for 
studying a variety of novel phenomena in quantum de- 
generate phases. In particular, s-wave Feshbach reso- 
nances have been widely used to study the crossover be- 
tween Bose-Einstein condensation (BEC) and Bardeen- 
Cooper-Schrieffer (BCS) superfluidity [I] |3 H G3 |B]. 
One of the greatest challenges in the field of trapped 
Fermi gases is the realization of superfluidity with pairs of 
atoms in a non-zero relative orbital angular momentum 
state [71 [HI IS]- In such systems, a rich variety of phases 
are expected to merge due to complex superfluid order 
parameters [TOl HU G2J Q2J [14] . Non- s-wave superfluidity 
has so far been realized in superfluid 3 He (p-wave) |15j. 
strontium ruthenate superconductors (p-wave) |16j. and 
high-T c copper-oxide superconductors (d-wave). 

Until now, 40 K and 6 Li have been utilized to study p- 
wave Feshbach resonances [HI [HI US] HH EH E21 [23] El] . 
The positions of p-wave Feshbach resonances of 6 Li have 
been accurately determined by Zhang et al. |17j and 
Schunck et al. [H]- Zhang et al. also detected the for- 
mation of p-wave Feshbach molecules of 6 Li in one of 
the three combinations of the two lowest atomic hy- 
perfine spin states, \F = 1/2, tuf — 1/2) (= |1)) and 
\F = 1/2, m F = -1/2) (= |2)), but molecules were 
not detected in the |1) — |1) state, which is the most 
promising candidate for a p-wave condensate because 
the dipolar loss is energetically suppressed. Recently, 
Fuchs et al. measured the binding energies and rela- 
tive molecular magnetic moments of |1) — |1), |1) — |2), 
and 1 2) — |2) molecules of 6 Li[TS]. Gaebler et al. suc- 
cessfully created p-wave Feshbach molecules of 40 K and 
observed the anisotropic energy release at the dissoica- 
tion process [20, 21J. The lifetime of 40 K molecules in 
both mi — ±1 and mi — states were found to be less 
than 10 ms, even on the bound side of the resonance due 
to dipolar relaxation into lower atomic spin states[2T]. 

While s-wave Feshbach molecules are highly stable 
against vibrational quenching near a Feshbach resonance, 
it is not clear how much p-wave molecules would suffer 



*Email address: muka@sogo.t.u-tokyo. ac.jp 



from inelastic collisions due to an intrinsically small in- 
teratomic separation caused by the centrifugal barrier. 
Therefore, studying collisional properties of p-wave Fesh- 
bach molecules is an important step toward realization of 
p-wave molecular condensates. In this Letter, we present 
elastic and inelastic collisional properties of 6 Li2 p-wave 
Feshbach molecules. To suppress the loss of molecules, we 
irradiated a resonant light pulse and blew away the resid- 
ual atoms from the trap immediately after the molecules 
were created [25l [26] [27]. The resultant pure molecular 
sample in the optical trap enables us to measure inelas- 
tic decay rates of molecules precisely. We measured the 
atom-dimer inelastic collision coefficients by keeping both 
atoms and molecules in the trap, and dimer-dimer inelas- 
tic collision coefficients by using a pure molecular sample. 
We also measured the dimer-dimer clastic collision rate 
from the thermalization rate of breathing-mode oscilla- 
tions of the cloud, which were induced spontaneously by 
the mismatch between the initial size of the atomic cloud 
and the equilibrium size of the molecular cloud [28 . 

Our experiments were performed using a quantum de- 
generate gas of 6 Li in the |1) and |2) spin states. In 
the present experiment, we evaporated the atoms in a 
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FIG. 1: (a) Time sequence of magnetic field to create 
molecules. Magnetic field is ramped adiabatically from above 
toward the Feshbach resonance (dashed line), (b, c) Absorp- 
tion images of dissociated molecules by the upward magnetic 
field ramp, (d, e) Atoms in a molecular gas detected at (d) 
or below (e) the resonance, (f) No atoms and molecules are 
detected when the magnetic field is off resonance. 
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FIG. 2: (a) Number of 1 1) — 1 1) molecules created as a func- 
tion of magnetic field detuning AB. The solid curve shows a 
Lorentzian fit to the data. Inset shows the time sequence of 
magnetic field for data taking, (b) Number of atoms remain- 
ing after a hold time of 25 ms as a function of AB. 



crossed beam optical dipole trap at 300 G, where the 
elastic collision cross section between |1) and |2) takes 
on a local maximum. After evaporation for 4 s, the mix- 
ture of 5 x 10 5 atoms in 1 1) and |2) with equal population 
and the temperature of T/Tp ~ 0.25 was prepared. The 
optical trap beams were orthogonal to each other on the 
horizontal plane: One, oriented along the x axis, has 
beam waist of 26 /im and the other, oriented along the 
y axis, has an ellipsoidal shape with a waist of 360 /im 
in the x direction and 170 /im in the z (vertical) direc- 
tion. The trap frequencies (u) x , uj y , lo z ) at the final trap 
depth are 2ir x (39, 760, 780) Hz. By stabilizing the cur- 
rents in the Feshbach coils, we suppressed magnetic-field 
fluctuations due to current noise down to 10 mG. 

To create |1) — |1) molecules, we first ramp the mag- 
netic field up to the Feshbach resonance of |2) — |2) at 
215 G to eliminate atoms in the |2) state. After prepar- 
ing the atoms in the pure |1) state, we create molecules 
by ramping the magnetic field to a value near the Fes- 
hbach resonance of the |1) — |1) state at 159 G. Due to 
the smallness of separation of Feshbach resonances for 
mi = ±1 and to/ = molecules [22 , we were unable to 
select a specific to; state. The time sequence used to 
create |1) — |1) molecules is shown in Fig. [lja). When 
the magnetic field reached the resonance, the sample was 
irradiated with a 25- /is pulse of resonant light ( "blast" ) 
to remove unpaired atoms from the trap. Molecules are 
relatively immune to the blast light due to the small- 
ness of the Franck-Condon factor for the optical transi- 
tion. After blasting the atoms, the magnetic field was 
ramped up above the resonance (Fig. [TJb,c)), kept at 
the resonance (Fig. pld)), or ramped down below the 
resonance (Fig. [TJe)j before taking absorption images. 
The blue, green, orange, and red curves in Fig. [lja) cor- 



respond to the magnetic-field sequence for the images in 
Fig. [TJb-e), respectively. Ramping the magnetic field to 
a value above the resonance as shown by the blue and 
green curves in Fig. [lja), allows us to clearly image the 
dissociated atoms(Fig. [TJb,c)). As shown by the or- 
ange and red curves in Fig. [TJa), the atoms in bound 
molecules are detected, albeit weakly, because of a finite 
Franck-Condon factor for the optical transition. When 
the magnetic field is kept away from the resonance by 
AB = +100 mG, no atoms were detected (Fig. [ljf ) ) , in- 
dicating the blast light completely removed the unpaired 
atoms. Therefore, the atoms imaged in Fig. [TJd,e) are 
not unpaired atoms but atoms that form molecules. As 
seen in Fig. [TJb-d), the higher the magnetic field was 
applied for dissociation, the more widely the cloud ex- 
panded during the same time of flight of 500 /is. This 
indicates that faster ramping of the magnetic field leads 
to higher dissociation energy. We produced |1) — |2) and 
1 2) — 1 2) molecules in a similar manner by using Feshbach 
resonances at 185 and 215 G, respectively. 

To investigate the range of magnetic field over which 
molecules can be formed efficiently, we ramped the mag- 
netic field from below up to the resonance. This avoided 
the adiabatic formation of molecules. As shown in the 
inset of Fig. |2^a), we first ramped the field across the 
resonance at a fast rate in order to minimize the atomic 
loss. We then swept the magnetic field in 1 ms to a 
value at which the number of molecules formed was mea- 
sured. After molecules were formed during a hold time 
of 500 /is, the blast light was turned on to remove un- 
paired atoms. Figure^ a) shows the number of molecules 
formed vs magnetic-field detuning AB, where AB is de- 
fined as positive when the magnetic field is higher than 
the resonance and the molecular state becomes quasi- 
bound. Throughout the paper, AB — is defined as 
the magnetic field at which the maximum number of 
molecules is created. At the peak, we can convert 15 
% of atoms to molecules, which is consistent with previ- 
ous experiments [T71 |2D]. Figure [2^b) shows the atomic 
loss measured after the 25-ms hold time for each mag- 
netic field. An asymmetry in the loss feature is clearly 
seen, as observed by other groups |T8| fT9] . Note that the 
bottom of the loss curve does not coincide with the peak 
of the molecular creation curve. The physical origin of 
this differecne is unclear and requires further study. 
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FIG. 3: Time evolution of T* . Thermalization time is esti- 
mated to be 20 ms. 
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In the process of creating molecules, breathing oscil- 
lations were excited spontaneously. This is because the 
molecules were created over the entire cloud of the Fermi- 
degenerate atoms, which is larger than the equilibrium 
cloud size of molecules [25]. Figure [3] shows the time 
evolution of T*(t), defined by T*(t) = muj 2 a 2 ms (t) /fee, 
where m is the mass of molecules, ui is 2ir times the trap 
frequency along the longest axis (x axis) of the trap, 
and <7 rms is the root-mean-square size of the molecular 
cloud. The size of the molecular cloud is largest when 
the molecules are created, and starts to shrink due to 
the optical trapping potential. The frequency of initial 
oscillations is evaluated to be 87(10) Hz, which agrees 
with twice the trap frequency of 39 Hz in the x direction. 
As seen in Fig. [3j T* (t) shows abrupt thermalizaion at 
t = 20 ms, by which we determined the thermalization 
time to be 20 ms. It is reported that it takes about 2.7 
collisions for the cloud to thermalize (29] [30] . Based on 
this, the dimer-dimcr clastic collision rate is estimated to 
be (20ms/2.7)~ 1 - 140 Hz. 

We measured the lifetimes of molecules in all three 
combinations of |1) and |2) states (Fig. E]). Figures Qa)- 
|4^c) show respectively the numbers of |1) — |1), |1) — |2), 
and 1 2) — |2) molecules remaining after the varying hold 
time. InQa), the number of molecule is greater than 
that in Fig. |2ja) due to different experimental conditions. 
Filled circles show the number of molecules vs hold time 
under the condition that only molecules are held in the 
trap, and open circles show the same quantity in the 
presence of unpaired atoms. In the latter case, the blast 
light was turned on after the hold time to ensure that the 
signal arises exclusively from molecules. The molecular 
lifetime is greatly suppressed in the presence of unpaired 
atoms. Magnetic-field detuning during the hold time is 
AB — — 175mG for all data points. 

The loss of molecules is analyzed using the following 
equation [27] 
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FIG. 4: Number of molecules vs hold time for |1) — |1) (a), 
1} — 1 2} (b), and |2) — |2) molecules (c). Filled circles show 
the decay data for molecules alone, and open circles show the 
decay data with both atoms and molecules. The solid and 
dotted curves are for pure molecules and for atom-molecule 
mixtures. 



where n a and n d are the density of atoms and that of 
molecules, respectively, a is the one-body decay rate of 
a molecule, and K ad and K dd are the loss coefficients for 
atom-dimer and dimer-dimer inelastic collisions, respec- 
tively. In fitting Eq. (JT|) to experimental data, we use 
the peak densities of atoms and molecules for n a and n d . 
We also assume that the spatial distribution of atoms 
obey an ideal Fermi-Dirac distribution in a harmonic 
trap, and n a is time-independent since the number of 
atoms far exceeds that of molecules [31]. For the profile 
of molecules, we assume a Gaussian spatial distribution 
rid = Ad/((27r) 3 / 2 (7 x a'j / cr z ), where cr x , cr y , and <r z are the 
root-mean square radii of the molecular cloud in the x, y, 
and z directions, respectively. Time independent cloud 
size is assumed in the analysis [27] [28] . 

The measured decay rates are summarized in Table [I] 
For |1) — |1) molecules with no unpaired atoms, the de- 
cay curve is well fitted by Eq. (JT|) without the one-body 



decay term. As mentioned above, the molecular cloud 
undergoes breathing-mode oscillations, and the molecu- 
lar density changes in time up to 20 ms. Therefore, we 
fit the decay data to extract Add only at t > 20 ms, 
until when the molecular cloud is expected to reach ther- 
mal equilibrium. For the |1) — |2) molecules alone in the 
trap, we used the whole data for the fit unlike the case 



TABLE I: One-body decay rate a, atom-dimer (Kad) and 
dimer-dimer (Add) inelastic collision coefficients for |1) — |1), 
1) - 1 2), and |2) - |2) molecules. 
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FIG. 5: Lifetime (open circles) and dimer-dimer inelastic 
collision coefficient (filled circles) for |1) — |1) molecules vs 
magnetic-field detuning. The dotted (dashed) lines show the 
averages of K^d (1/e lifetime) for AB < — 50mG. 



of |1) — |1) molecules because of the short lifetime. Since 
|2) is not the lowest energy state, |1) — |2) molecules are 
expected to undergo dipolar decay into the lower spin 
state. Therefore, both a and K^d contribute to the de- 
cay. With the initial density rid = 4.8 x 10 11 cm~ 3 , the 
decay rate Kdd^d is evaluated to be 390 Hz, which is 15 
times larger than the one-body decay rate a. For |2) — 12) 
molecules alone in the trap, an exponential decay fits the 
data well, indicating that the one-body decay is domi- 
nant. While we are able to extract K a d for the |1) — |2) 
molecules from the decay data, it is difficult to extract 
Kad for 1 2) — 1 2) molecules due to the low signal-to-noise 
ratio. Here, n a is estimated to be 4.5, 8.0, 4.5 x 10 12 cm -3 
for Fig. |3Ja)-|4jc) , respectively, nd is estimated to be 0.9, 
4.8, 1.1 xlO"cm- 3 for |1) - |1) (at t = 20 ms), |1) - |2) 
(at t = ms), |2) — |2) (at t = ms), respectively. 

We measured the AB dependence of the 1/e lifetime 
and -Kdd of the |1) — |1) molecules. The 1/e lifetime 



of molecules is plotted in Fig. [5] (open circles). For 
— 340mG < AB < — 50mG, the lifetime is independent 
of AB and similar to the thermalization time of 20 ms. 
The filled circles in Fig. [5] show Kdd evaluated from the 
decay curve in the same manner as in Fig. Qa). Since 
the lifetime is expected to become shorter when the mag- 
netic field crosses the resonance, our uncertainty in the 
resonance position could be as large as 50 mG. We eval- 
uate and plot Kdd in Fig- [5] only in the range between 
—340 mG and —50 mG because neglecting the one-body 
decay term is valid only in the bound-molecular side of 
the Feshbach resonance. 

In conclusion, we studied the elastic and inelastic col- 
lisional properties of 6 Li 2 p-wave Feshbach molecules. 
From the thermalization time of the breathing-mode os- 
cillation, we evaluated the dimer-dimer elastic collision 
rate. From the loss measurements expressed as a function 
of the hold time, the dimer-dimer inelastic loss coefficient 
for |1) — |1) and |1) — |2) molecules and the one-body de- 
cay coefficients for |1) — |2) and |2) — |2) molecules were 
obtained. We also evaluated the atom-dimcr inelastic 
collision rates for |1) — |1) and |1) — |2) molecules. In 
the present experiment, we created 7 x 10 3 molecules at 
a temperature of 1 /iK after thermalization. Provided 
molecules are distributed equally among three different 
mi states, the phase space density of the molecular gas 
is estimated to be 4 x 10~ 3 . Our measurements show 
that the ratio of elastic to inelastic collision rate is about 
five, indicating that a conventional evaporative cooling 
would not be efficient [32] . To overcome this, one could 
use Feshbach resonance to enhance the dimer-dimer clas- 
tic collisions. One can also utilize the optical lattice to 
single out a specific mj state [H] which should have a 
higher ratio of elastic to inelastic collision rate. 
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